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FROST FORMATION MODELING APPLIED TO VARIOUS DATA SETS 

SECTION 1 
INTRODUCTION 

Frost formation commonly occurs when a surface in the 
presence of humid air is cooled to temperatures below 0*’C and 
below the dew-point. In the field of refrigeration and air- 
conditioning this phenomenon has an adverse effect upon heat 
transfer and pressure loss^ thus leading to lower efficiencies. 

In aerospace technology frost creates safety hazards. For 
instance, in the northern climates, aircraft takeoffs can be 
hazardous if nocturnal frost is present on the airfoil. Another 
technology concerned with frost is that of the space shuttle. A 
frost layer on the external cryogenic tanks of the space shuttle 
could shear off at launch and damage the orbiter tiles. The 
emerging cryogenic technology in other industries is very likely 
to experience similar concerns with frost formation. 

Frost formation is a complicated transient process in which 
a variety of heat and mass transfer mechanisms are at work 
simultaneously. It can be characterized as having three regimes. 
In the first regime, the initial frost layer can begin in one of 
two ways. In the first case initial condensation occurs at 
nucleation sites on the wall resulting from a critical super- 
saturation. In the other case (for a very cold wall) boundary 
layer fogging occurs and the fog becomes the major source for 
water droplet condensation on the wall. Thus the "frost-point" 
temperature (i.e., the temperature at which frost actually begins 
to form) lies somewhere between the fog-point temperature and the 
dew-point temperature, depending on nucleation sites on the wall 
and on the wall temperature. The "frost-point" temperature is an 
especially important consideration for nocturnal frost formation. 
Some initial frost properties such as density and thickness need 
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to be calculated, because of the processes described above, in 
order to move on to the next regime of frost formation. 

The second regime is dominated by vapor diffusion from the 
air boundary layer into the macroscopic frost layer. This pro- 
cess increases the average frost density. As the frost grows in 
thickness, the insulating effect of the frost layer causes an 
increase in the frost surface temperature. This in turn causes a 
decrease in the mass flux, which consists of both vapor flux and 
fog droplets flux in the boundary layer, as well as a decrease in 
the heat flux. The frost layer thickness also increases when the 
mass flux from the boundary layer is larger than the diffusional 
mass flux entering the frost layer. Boundary layer fogging may 
disappear (depending on the frost surface temperature) since the 
fog droplets are assumed to freeze on the frost surface. 

Generally, since the thermal conductivity of the frost increases 
along with frost density, this regime is characterized by a feed- 
back effect in which the frost surface temperature slowly 
increases as the frost layer grows. 

After a while a third regime is reached. The frost surface 
temperature becomes a constant, thus giving a steady-state heat 
flux in the boundary layer. The phenomena of liquid water 
seepage, in addition to the water vapor diffusion, from the frost 
surface into the interior of the frost need to be considered in 
order to explain the constant frost surface temperature. An 
effect that remains mysterious, but has been acknowledged by many 
investigators is that, at least on a flat plate, the frost den- 
sity at any time tends to be uniform throughout the frost layer, 
although the density increases with time. The spatially uniform 
frost density can not be completely explained by a vapor dif- 
fusion model. One possible explanation assumes a thermal dif- 
fusion process of water or ice particles within the frost layer. 

No one has yet constructed a model based on such a process, 
however. The assumption of a spatially uniform frost density has 
been utilized in the frost formation model described in this report. 
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In classic work, Brian*, et al. [1] and Brian, et al. [2] 
chose a diffusional model approach to calculate the frost density 
and thickness as a function of time. They used a heat balance at 
the frost surface to calculate the frost surface temperature. 

Then by assuming the frost surface is saturated and the frost 
density is spatially uniform, the rate of frost density change 
was related to the temperature gradient at the frost surface. 
Finally, the difference between the boundary layer vapor mass 
flux and the diffusional mass flux into the frost could be used 
to calculate the frost thickness. 

When the Brian model is closely examined, one discovers 
certain limitations. First, in the heat balance equation, the 
frost thermal conductivity is only applicable at wall tem- 
peratures of 80®K and with a high initial saturation ratio, as 
demonstrated by Dietenberger [3]. Second, the heat and mass 
transfer coefficients of the boundary layer are not very well 
known, particularly in relationship to the roughness of the frost 
surface. Chen and Rohsenow [4] showed that the heat and mass 
transfer coefficients are affected by the frost surface rough- 
ness. Third, the calculation of frost density and thickness is 
dependent on an effective numerical integration scheme which also 
requires fairly good initial values of frost density and 
thickness. Finally, there are no provisions for the quasi-steady 
state phase, where the frost surface temperature reaches a 
constant. Because the diffusional approach lends itself to time- 
varying meteorological conditions, Jones and Parker [5] applied 
Brian's et al . method and achieved partial success in predicting 
frost thickness trends, even when humidity and/or velocity 
changed during an experimental run. 

Due to the complicated nature of frost growth, other early 
investigators attempted an empirical approach to the frost den- 
sity prediction without relying on a physical model. A signifi- 
cant effort was made by Biguria and Wenzel [6] to correlate frost 
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density on a flat plate with air stream velocity, air stream 
humidity, plate temperature, nature of the boundary layer, and 
the frost surface temperature. Lotz [7] used a different data 
set to derive yet a different correlation of the frost density 
with frost thickness, air stream velocity, air stream temperature, 
and the vapor pressure difference between the air and frost. 
Lastly, Parish and Sepsy [8] used a simple correlation of the 
frost density with frost surface temperature to calculate the 
frost growth on a cylinder in conjunction with a laminar boundary 
layer/potential flow method. Without any explanation of a physi- 
cal model of frost density, it was not imperative to apply the 
frost density correlations to the model of frost growth in 
general . 

Recently, some investigators who used the quasi-steady 
state assumption were able to postulate a physical model of frost 
thickness growth. In the quasi-steady state assumption the frost 
surface temperature or the heat flux eventually reaches a 
constant value, although the frost density and thickness will 
continue to increase. If the frost grows long enough, the frost 
surface eventually reaches 0®C, providing that the air stream 
temperature and the dew-point are above freezing temperature. In 
experiments for forced convection over a cylinder, Schneider [9] 
observed that frost thickness tends to be independent of the 
Reynold's number and the vapor pressure difference between the 
air stream and the frost. Consequently, he formulated frost 
thickness according to the principles of crystal growth. Treating 
a frost crystal as a vertical cylinder, the heat flux of conden- 
sation was equated with the heat flux of conduction through the 
cylinder. This resulted in a frost thickness growth written as 
the square roots of the time and of the temperature difference 
between the frost surface and the wall. This basic model was 
empirically corrected by including the effects of supersaturation 
and the condensation soaking of the porous frost structure. 
Schneider then used five other references containing experimental 
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data on cylinders to claim that the generality of his correlation 
for the frost thickness was independent of the geometrical shape 
of the flow channel as well as of the kind of flow. In this 
paper we will show that Schneider's correlation is valid for his 
own data set, but not in agreement with other data selected for 
review. 

In a more basic development. White [10] used thermal con- 
ductivity as a linear function of the frost density and as an 
arbitrary function of the frost interior temperature. By 
invoking the quasi-steady heat flux assumption and integrating 
the frost thermal conductivity across the frost thickness, the 
overall frost density became a linear function of the frost 
thickness. This operation was then combined with the mass 
balance at the frost surface, and spatially invariant frost den- 
sity was assumed, resulting in frost thickness as an analytic 
function of time, mass flux to the frost surface, the initial 
frost density and thickness, the wall and frost surface tem- 
peratures, and the coefficients of the thermal conductivity. 

White found that for the range of experimental data observed, 
there was a very small dependence on the wall and frost surface 
temperatures and the coefficients of the thermal conductivity. 
From these observations he concluded that half the mass flux is 
deposited on the frost surface and the other half goes into the 
frost to increase the frost density. We note that these correla- 
tions were limited to low frost densities and are only applicable 
to the quasi-steady regimes of frost growth. 

This paper describes a numerical model that has been used 
to calculate the frost properties for all regimes of frost 
growth. In the first regime of frost growth, the initial frost 
density and thickness was modeled from the theories of crystal 
growth, similar to Schneider's approach. The "frost-point" tem- 
perature was modeled as a linear interpolation between the dew- 
point temperature and the fog-point tem.perature , based upon the 
nucleating capability of the particular condensing surface. 
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For a second regime of frost growth, the diffusional model 
of Brian, et al, was adopted with the following enhancements: 
the generalized correlation of the water frost thermal conduc- 
tivity from Dietenberger was applied to practically all water 
frost layers in place of Brian's, et al. correlation, being care- 
ful to ensure that the calculated heat and mass transfer coef- 
ficients agreed with experimental measurements of the same 
coefficients. In particular, the heat transfer coefficient was 
enhanced due to the roughness of the frost layer and the mass 
transfer coefficient was enhanced due to the boundary layer 
fogging. Finally, an efficient and accurate numerical integra- 
tion scheme was developed to calculate the frost surface tem- 
perature and the frost density and thickness as a function of 
time. The numerical scheme responded to changes in environmental 
conditions, with good agreement to Jones and Parker's data. 

In the third regime of frost growth, the frost surface 
temperature was held fairly constant by a procedure described as 
the water seepage model, which is explained later in this paper. 
Agreement with experimental data for the quasi-steady regime was 
superior to that obtained using Schneider's or White's 
correlations. 


6 



SECTION 2 

THE FROST INITIATION MODEL 


The wall temperature at which the frost begins to form is 
the frost-point temperature/ which is determined by the equation/ 

% = ’"dp 

The dew-point temperature and the fog-point temperature 
equations are presented in Appendix I. 

Once the frost begins to form/ we assume the initial frost 
layer will grow according to the principles of crystalization. 
Following Schneider [9] the frost layer is modeled as cylindrical 
ice needles growing away from the wall. The one-dimensional 
conductive heat flux through the ice needle is assumed balanced 
by the various heat flux contributions from the air to the ice 
crystal. The surface temperature of the ice needle is assumed to 
be at the dew-point temperature but not greater than the freezing 
temperature. The heat flux terms are different than those assumed 
by Schneider. The steady state conductive heat flux through the 
ice crystal is given by 


Q 


o 


630 ^ 
T dx 


630 


Jin(T 


/T ) 


(2) 


This heat flux is balanced by the heat flux from the air to the 
ice as/ 


Q “ ^ ^ p Ti 

o o s^lV 


u) - 0 ) \n /dx 
®a 


a Wi 


H 


V -r 

(3) 


(4) 


In Schneider's approach, n is 1/2 and is zero, which 
worked well for thick frost layers. The classical nucleation 
theory requires n=l for an ice needle isolated from other 
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needles r as would be expected in an initial frost layer. As will 

be demonstrated later , H turns out to be small compared to Q . 

o o 

Solving for the ice needle height or the initial frost thickness, 
Xs^f we obtain by integration from equations (2) and (3) and with 
n=l. 



0^3“!) 



t 



(5) 


where 


= 630tn(Tg /Tj 

o 


( 6 ) 


During the short time needed for the initial frost layer 
to develop, the mass flux to the frost layer from the air is 
approximated as, 


h 


m 


( u -u» ) 
V a wy 


(l-B)p 



(7) 


Combining equations (5), (6), and (7), we obtain the frost den- 
sity as a function of x , as, 

s_ 


Pf 


o 




/C h L 
/ o m s 

V 




In order to estimate a value for x so that can be 

s ^f 

calculated, we use the approximation 

1 ^ 

Jin(l-x) « "55 “2^ 

in equation (5) to explicitly solve for x as 


( 8 ) 


X 


s. 



Zn(Tg^/T^) 


V2 


(9) 


8 



and 



1260 t 


Vs"! 


( 10 ) 


The terms C and C, are fitting constants obtained from the data 
o 1 

reviewed in this paper. The nominal values for and are 
15.56 and 0.001. Equations (8) and (9) allow us to solve for the 
initial frost density and thickness; that is, the thickness and 
density at the end of the crystal growth dominated regime. These 
values are inputs to the diffusion dominated regime modeling. It 
turned out that the values of initial frost density and thickness 
did not have to be very accurate because the diffusion growth 
model generates most of the frost layer. Thus at times far 
removed from the initial time, the frost thickness and density 
are not very sensitive to the intial values of frost density and 
thickness . 
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SECTION 3 

THE FROST FORMATION MODEL 


The basic approach for modeling frost growth in the 
diffusion dominated regime follows that of Brian, et al. except 
for differences in the various terms of the equations. When the 
heat flux, which includes latent, convective, and radiative 
transport from the air to the frost surface, is made equal to the 
conductive heat flux through the frost layer at the frost surface, 
the temperature of the frost surface can be calculated. The 
frost surface temperature will change with time because the frost 
layer becomes denser and thicker, which in turn affects the com- 
putations of the heat flux in the frost layer. During the pro- 
cess of frost formation it is assumed that the part of the water 
vapor and all of boundary layer fog droplets transported to the 
surface freeze at the surface, thus increasing the thickness of 
the frost. The remaining water vapor is diffused into the 
existing frost layer before it freezes. Since the frost density 
is assumed to be spatially invariant in a direction normal to the 
plate, the water vapor diffusion flux entering the frost surface 
from the surrounding air is given- by 


^ 3 (1-B) 
ds 3t *^I^s 


( 10 ) 


where m, , p_ , and x stand for water vapor diffusion mass flux 
at the surface, ice density, and frost thickness respectively. 

B is the porosity of the frost. Also, the water vapor mass flux, 
m^ , will obey the diffusion equation and will be driven by the 
temperature gradient throughout the frost layer. 

The water diffusion flux is given by equation (8) from 
Dietenberger , 

P 

Ji - DB V 

"‘d = (l-x)Xs ^^2 

V 
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dT 

gradient, 


Equation (11) can be evaluated at the frost surface if the 

the frost surface temperature 
, and the frost thickness x are known. Combining 
) and (11) allows the frost density to be solved for 


frost surface temperature, , 


dx 

equations (10' 


as a function of time. 


To obtain a value for frost surface temperature, T, 


the 


quasi-steady-state heat equation for the frost layer is to be 
solved. That is, the equation 


trdT 

^dx ” ^ 


( 12 ) 


with the frost thermal conductivity, K, given by equation (23) 
from Dietenberger (reproduced in Appendix II) is to be solved 
with the boundary conditions 


where 


T = T at X = x„ and T =T at x = 0, 
s s w 


q = h„ (T -T ) + h L (uj - lo ) 
^ H as msa s 


4 4 

+ ea (T^^ - T^^) 
a s 


(13) 


This gives 


/ K dT = x q 

T ® 

w 


(14) 


The temperature gradient at the frost surface is given by 


dT ^ q 

dx s K(T ) 
s 


(15) 


The thickness of the frost can be computed from the ice 
portion of the frost (l-B)p x which is directly related to the 
mass transfer coefficient and the water vapor humidity by 


d[(l-B) pjXg] 
dt 


h 


m 


(ai^-cos) 


( 16) 
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since these equations cannot be solved analytically, numerical 
techniques are used to solve these equations to obtain the 
thickness and the density of the frost. 

The numerical scheme begins by assuming a small initial 
value of frost density, pfQ and frost thickness, XgQ, obtained 
from equations (8) and (9). 

The frost surface temperature, Tg, is obtained from the 
monotonic nonlinear equation (14). In order to solve equation 
(14) in an accurate and efficient manner, the frost thermal 
conductivity, K, in equation (14) is first interpolated by a 
parabolic function in a temperature range AT and integrated to 
get the first term in equation (14) (i.e., AT is often greater 
than Tg Then Tg in equation (14) is solved for by an 

interval-halving-iterative technique so that the left term of 
equation (14) is equal to the right terms. In this way suc- 
cessive evaluations of the complicated thermal conductivity 
expression, K, and numerical integrations can be completely 
avoided when iterating to get Tg. The temperature gradient at 
the frost surface is obtained by equation (15). The value for 
the frost density is predicted for a later time, t^+i^ by com- 
bining equations (10) and (11) to get 3pf/3t, which is solved by 
using a second order Runge Kutta "initial" scheme and by a second 
order predictor-corrector "takeover" scheme. An accurate predic- 
tor equation was adapted from White's [10] analytic form for the 
frost density. Finally, the value for the frost thickness is 
predicted for the same time as the frost density (i.e., at t^+i) 
by integrating equation (16), where the rectangular rule is used 
for the integration, in acknowledgement of the fact that hj^ and 
(Og are slowly changing functions of time. At the same time, a 
new frost surface temperature is calculated by returning to 
equation (14), incrementing n, and using the new values of frost 
density, pf, and thickness, Xg. 

The regime of the quasi-steady state where the frost sur- 
face temperature is held constant is modeled in the following 
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way. If at a particular time, t^, the frost surface temperature, 

T , as solved by the previously described method, goes above a 
s 

critical temperature, hereafter known as the water seepage tem- 
perature, then the numerical scheme is altered to let the frost 

thickness be a constant in the time interval t to t , and the 

n n+1 

frost density is computed by the integration of equation (16) 

over the interval t^ to Thus liquid water seeps into the 

frost layer, increasing its density but not affecting the 

thickness. The time is incremented and the frost surface 

temperature, T , is solved from equation (14) and is again com- 
s 

pared to the water seepage temperature, and so on. Altering the 
frost growth in this manner has the effect of keeping the frost 
surface temperature constant. Generally the water seepage tem- 
perature tends to be lower than or equal to 273.16*K in order to 
gain agreement with the data presented later. It is postulated 
that some isolated regions of the frost surface are at 273.16®K 
while the other regions are below the water seepage temperature, 
which when the region temperatures are averaged out gives the 
frost surface temperature. Thus, for very cold wall temperatures 
the water seepage is expected to be somewhat less than freezing 
while for wall temperatures approaching freezing the water 
seepage temperature would nearly equal the freezing temperature. 

The frost formation model as described above need not be 
restricted to flat plate frost formation. The model can be 
applied to any geometrical surface including curved surfaces 
as long as the wall temperature and the heat and mass transfer 
coefficient distribution can be provided. 
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SECTION 4 

COMPARISON TO EXPERIMENTAL DATA 

The validity and generality of the frost formation model 
has been demonstrated by comparison with several different data 
sets appearing in the literature. Each set of published data 
required separate derivations of the convective heat and mass 
transfer coefficients depending on whether the air flow was 
forced or natural, in a duct, on a cylinder, laminar or 
turbulent, etc. Each data set contained the appropriate values 
needed for input into the frost formation model. These values 
include wall temperature, air temperature, relative humidity, 
airspeed, and the apparatus dimensions. In addition, many 
published data sets also presented measurements of both frost 
weight and thickness. Other published data just contained the 
frost thickness. For the latter data sets, a very accurate mass 
transfer coefficient was needed to compensate for the lack of a 
frost weight measurement. 

For each data set presented, the analysis of the convective 
heat and mass transfer coefficient will be discussed first. Then 
a table of input values to the frost formation model will be pre- 
sented for that data set. Finally the corresponding frost model 
prediction of frost weight and thickness will be compared to the 
experimental observations and an evaluation provided. The 
published data examined were from Brian et al, [2], Brian, et al , 
[1] , Yamakawa, et al. [11] , Jones and Parker [5] , Okino and 
Tajima [12], Schneider [9], Andrichak [13] and White [10], 

4.1 TURBULENT FORCED CONVECTION IN A DUCT 

Most of the data sets analyzed were frost formation experi- 
ments performed in a duct with turbulent forced convection. 

These included Brian et al, [2], Brian, et al. [1], Yamakawa, et 
al . , and Jones and Parker, 
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In Brian r et al. [2] and Brian et al. [1] experimental 
runs, the boundary layer temperature profiles were not fully 
developed but were fully turbulent. In other words the entrance 
effects of the duct were affecting the heat and mass transfer but 
were not affecting the momentum transfer. For the Yamakawa, et 
al. experiments the boundary layer temperature profiles were both 
fully developed and fully turbulent. The air velocities and duct 
dimensions of the Jones and Parker experimental runs are similar 
to Yamakawa et al., thus leading to similar conclusions as 
Yamakawa et al., concerning boundary layer temperature and velo- 
city profiles. The simplicity of the flow conditions allows one 
to use the Nusselt and Sherwood number correlations commonly used 
in the general literature for the calculation of heat and mass 
transfer coefficients. In particular, Yamakawa et al. used the 
popular Colburn equation, 

Nu = 0.023 Re°*®pr^/^ = h„D /k . (17) 

n 6 cl 

for the Nusselt number, Nu, as a function of Reynold's number. 

Re, based on the hydraulic diameter, D^, and the Prandtl's 
number, Pr. This curve is plotted in Figure 1. 

The Sherwood number is usually obtained via the 
Chilton-Colburn analogy, from the Nusselt number as 

W-? ^ ° 

Sh = Nu(Sc/Pr)-^' (18) 

*^a a 

or equivalently equations (17) and (18) can be rewritten in terms 
of heat and mass transfer coefficient ratio as 

hfl/hm = Cp(Pr/Sc)"^/^ . (19) 

These relationships assume that the duct walls are smooth 
and clean and that the boundary layer temperature and velocity 
profiles are fully developed and fully turbulent. In the 
Yamakawa, et al. case these assumptions were questionable and 
therefore Yamakawa et al. checked experimentally to make 


15 




Correlation of the Experimental Nusselt's Number Containing 
the Heat Transfer Coefficient, h*, versus Reynolds Number. 


Figure 1. 


sure that Colburn's equation for the Nusselt number equation 
was still valid. Their results in Figure 1 show that the right 
hand term of equation (17) must be increased by a factor 1.95, 
which then must be explained. One possible explanation is an 
increase in heat transfer due to the surface roughness of frost. 
It was determined that two other possible causes of heat transfer 
enhancement — the boundary layer fogging effect on temperature 
profile, and the boundary layer transpiration as related to the 
vapor mass flux divided by the bulk flux — are nil. 

There is evidence that the heat transfer coefficient is 
affected by the frost surface roughness. This problem of frost 
roughness was taken up by Chen and Rohsenow [4] who did frost 
formation in cylindrical tubes. They attempted to calculate the 
aerodynamic roughness height of the frost layer and then calcu- 
late the corresponding increase in the heat transfer coefficient. 
Using the relative correction to hg from Yamakawa as 1.95, one 
can use the equations from Chen and Rohsenow to show that the 
approximate aerodynamic roughness height of Yamakawa 's data can 
be related as kg « Xg/0.6. This is the value of aerodynamic 
roughness due to frost recommended by Langston [15] in his stu- 
dies of hoar frost on airplane wings. Thus surface roughness 
appears to be a reasonable explanation of enhanced heat transfer 
coefficients . 

Yamakawa, et al . also experimentally checked to see if 
the Chilton-Colburn ' s analogy between the heat and mass coef- 
ficients was valid for their experiments. Figure 2 shows their 
plot of the experimental heat and mass transfer coefficients. 

They did not relate their data to relative humidity or Reynolds 
number as was done in Figure 1, Yamakawa 's data points lie above 
the relationship given by equation (19). On the other hand, the 
rough wall heat transfer coefficient is greater, by a factor of 



Kcal/ni 


h* 9/m^ • sec 
m 


12.5 25 37.5 50 62.5 75 



h* kg/tn^ • hr 


Figure 2, Correlation of the Experimental Heat Transfer Coefficient, h , 
versus the Experimental Mass Transfer Coefficient, h*. ^ 

(Data obtained from Yamakawa, et al.). ^ 



. w/n)^ 



analogy relationship to the rough wall heat transfer coefficient. 
This resulted in the upper line in Figure 2. The lower line was 
obtained by using the constant 1.5 instead of 1.95. 

If one proposed an enhancement in the mass transfer coef- 
ficient due to the boundary layer fogging as an explanation for 
the lower line in Figure 2, then the mass transfer enhancement 
averages out to 1.95/1.5 = 1.3. This is approximately the value 
typically predicted when boundary layer fogging is considered, as 
is done in the model (see Appendix I). 

Thus we conclude that to make sense of Yamakawa's data for 
the heat and mass transfer coefficients, the smooth wall heat 
transfer coefficient must be increased by the roughness factor 
1.95 and the smooth wall mass transfer coefficient increased by 
the mass transfer enhancement equation due to the boundary layer 
fogging. From Street [16] and Yaglom and Kader [17] we find 
that, at least for the fully rough flow, a basic analogy exists 
between the heat and mass transfer coefficients. This means if 
the Yamakawa et al. data is in the fully rough regime, then our 

I 

conclusion could be invalid. However, using equations for the 

friction coefficient from Chen and Rohsenow and a roughness 

height of k = X /0.6, we find that the Yamakawa et al. data was 
s s 

clearly in the transitionally rough regime, where a heat and mass 
transfer analogy does not apply. 

A similar analysis of the heat and mass transfer coef- 
ficients for Brian et al. [2] data shows the same type of 
conclusion, with one additional factor necessary to correct for 
the thermal entrance length of the duct. The single solid square 
point in Figure 1 replotted from Brian et al . [2] is the smooth 

wall Nusselt number measurement, which is 2.24 higher than the 
Colburn equation at the same Reynold's number. This is due to 
the thermal entrance effects of the duct as discussed by Brian et 
al. [2]. Also non-dimensional i zed and replotted from Brian, et 
al. [2] is the analogy of the Nusselt number via equation (18), 
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the Sherwood number times the one-third power of the Lewis 
number. These are given by the solid dots in Figure 1, which 
then shows that the mass transfer coefficient via the Sherwood 
number corresponds to the smooth wall case, allowing for the mass 
transfer enhancements due to boundary layer fogging. 

On the other hand, the heat transfer coefficient via the 
Nusselt number (corresponding to the solid triangle in Figure 1) 
shows a correlation of 3.36 higher than the Colburn's equation 

4 

(ignoring the Reynolds number at around 1.5 x 10 ) or 1.51 
greater than the smooth wall^ Nusselt number correlation. Using 
the factor 1.51 for the relative increase in the heat transfer 
coefficient due to frost roughness and using equations from Chen 
and Rohsenow we also find that the Brian et al . data are clearly 
in the transitionally rough regime. Thus these data sets, in the 
transitionally rough regime, demonstrate the following behavior of 
the heat and mass transfer coefficients; the mass transfer coef- 
ficient shows a smooth wall behavior in the transitionally rough 
regime and the heat transfer coefficient shows a rough wall beha- 
vior in the same regime. 

At this point, the only other information needed for input 
to the numerical frost formation model of turbulent forced con- 
vection in a duct is the specification of ambient conditions such 
as Reynold's number, absolute humidity, wall temperature, and air 
temperature. Since the frost formation model was adapted from 
Brian, et al. [1] we first compared the numerical predictions 
with Brian, et al. [2] data. Initially their time frames were 
translated to correspond to the zero time of frost formation 
rather than the zero time of their run. Table I provides a list 
of input data into the numerical frost formation model and 
Figures 3 and 4 show the corresponding calculations of frost 
thickness and the frost mass per area (here defined as frost 
weight) respecti\'ely from zero time of frost formation. The 
frost thickness data in Figure 3 shows a lot of scatter but on 
the average the model predicted the frost thickness fairly well. 
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TABLE I 

DATA INPUT TO THE FROST FORMATION MODEL FOR COMPARISON 
WITH BRIAN ET AL. [2] DATA FOR FORCED CONVECTION IN A DUCT 


Computer 
Curve No. 

Experimental 
Graph Symbol 

Reynolds 

Number 

Absolute 

Humidity 

Wall 

Temperature 

(°K) 

Air 

Temperature 

(°K) 

1 

o 

5100 

0.0031 

79.83 

307.05 

2 

A 

5790 

0.0031 

79.83 

297.60 

3 

D 

8000 

0.00357 

79.83 

297.60 

4 

• 

8430 

0.0035 

79.83 

284.27 


0.9 


to 

to 



Time (Minutes) 

Figure 3. Frost Thickness versus Time for 
(Reference 2) . 




Weight (q/cm^) 







Then in Figure 4, we see that the frost weight is predicted very 
well. The Brian, et al, [1] data reported more measurements of 
frost thickness and weight. This time the smooth wall heat and 
mass transfer coefficients were increased by 1.27 rather than 
2.24 due to the thermal entrance effects. The roughness effects 
increased the smooth wall heat transfer by 1.76 rather than 1.51 
as in Brian, et al. [2]. The input data are shown in Table II 
and the corresponding frost thickness and weight comparisons are 
shown in Figures 5 and 6 respectively. There is less scatter in 
the frost thickness measurements and a very good prediction of 
the data. There is a similarly good result for the frost weight 
in Figure 6. 

For the Yamakawa et al . reported measurements of frost 
thickness and weight, the input data are shown in Table III and 
the corresponding frost thickness and weight are shown in Figures 
7 and 8 respectively. The curves 1, 2, and 4 in Figure 7 show 
good agreement with the frost thickness data, while curves 3 and 
5 show the result of applying the water seepage model at the 
water seepage temperatures of 268. 9“K. The water seepage model 
was applied in retrospect after evaluating other published data. 

A rule of thumb of determining when the water seepage model needs 
to be applied is in the observation of when the frost thickness 
data suddenly deviates from the basic predicted curve. A more 
direct method is to use an experimenter ' s own observations of 
cycling frost surface temperature, thus indicating water seepage 
effects. As more data is analyzed, the use of the water seepage 
model will become more objective. Lastly, Figure 8 shows good 
agreement of predictions with the data of the frost weight as a 
function of time. 

The last example of frost formation in a duct with tur- 
bulent forced convection is from Jones and Parker. They only 
measured the frost thickness. No measurement was made of the 
frost weight or the heat and mass transfer coefficients. 

However, their apparatus and ambient conditions were similar 
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TABLE II 

DATA INPUT TO THE FROST FORMATION MODEL FOR COMPARISON 
WITH BRIAN ET AL. [1] DATA FOR FORCED CONVECTION IN A DUCT 


Computer 
Curve No. 

Experimental 
Graph Symbol 

Reynolds 

Number 

Relative 

Humidity 

Wall 

Temperature 

(»K) 

Air 

Temperature 

(“K) 

1 

0 

14500 

26.1 

79.83 

297 

2 

A 

9316 

26. 2 

79.83 

297 

3 

D 

5603 

26.0 

79.83 

297 

4 

• 

14750 

20. 8 

79.83 

297 

5 

A 

9247 

15.8 

79.83 

297 

6 

- 

5625 

16.2 

79.83 

297 






Figure 6, Frost Weight versus Time 
(Reference 1). 
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TABLE III 

DATA INPUT TO THE FROST FORMATION MODEL FOR COMPARISON 
WITH YAMAKAWA ET AL. [10] DATA FOR FORCED CONVECTION IN A DUCT 


Computer 
Curve No. 

Experimental 
Graph Symbol 

Air 

Velocity 

(m/s) 

Relative 

Humidity 

Wall 

Temperature 

(°K) 

Air 

Temperature 

(°K) 

1 

□ 

1.6 

56.2 

251 

284 

2 

A 

5.0 

56.2 

251 

284 

3 

o 

5. 0 

70.0 

251 

284 

4 

A 

7.3 

56.2 

251 

284 

5 

• 

7. 3 

97.4 

251 

284 

Note ; To 

calculate Reynolds 

Number, the 

hydraulic 

diameter is taken 

as 


D^ = 20(1.7 - x^)/(11.7 - x^) cm. 

6 S S 




Figure 7. Frost Thickness versus 
(Reference 10). 



Seepage begins =jsl- 




Time (minutes) 


Figure 8, Frost Weight versus Time for Yamakawa, et al. Data 
(Reference 10). 


enough to Yamakawa et al. that the correlations for the Nusselt 
and Sherwood numbers were extended to Jones and Parker's data as 
a best guess. Indeed ^ the Jones and Parker data warranted some 
attention because the ambient conditions were time-dependent, 
unlike other published data. The time-dependent ambient con- 
ditions are indicated in Table IV and the corresponding frost 
thicknesses in Figure 9. The very good comparison with the frost 
thickness data for the first sixty minutes prior to any ambient 
condition changes or any water seepage effects confirm the heat 
and mass transfer coefficient calculations. Curve 1 in Figure 9 
shows a discontinuous change at 60 minutes as it should and con- 
tinues to follow the data without the inclusion of water seepage 
effects. Curve 2 in Figure 9 also shows a dincontinuous change 
at 60 minutes and follows the data until 90 minutes without the 
use of the water seepage model. When a water seepage temperature 
of 268. 7"K is used, good agreement is obtained from 75 minutes to 
120 minutes. Curve 3, with a water seepage temperature of 
272. 0®K, is very interesting in that it levels out after 105 
minutes just as the data does. Curve 4, with water seepage tem- 
perature of 269. 3® K, shows a similar pattern to Curve 2 in pre- 
dicting the data. The dashed lines show the effect of setting 
the water seepage temperature at the maximum possible value of 
273.16“K. Thus the results are sensitive to selection of the 
water seepage temperature. 

4.2 LAMINAR AND TURBULENT NATURAL CONVECTION ON A VERTICAL 

PLATE 

The results from Nakamura [18] show that in laminar natural 
convection on vertical plates, the surface roughness has no 
effect on the heat transfer coefficient because the sensible heat 
flux from the air to the frost is mostly conductive, rather than 
convective. However, for natural convection, the momentum, heat 
and mass transfer mechanisms are coupled together in the boundary 
layer. The analysis for this problem is taken from Okino and 
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TABLE IV 

DATA INPUT TO THE FROST FORMATION MODEL FOR COMPARISON 
WITH JONES AND PARKER DATA [4] FOR FORCED CONVECTION IN A DUCT 


Computer 
Curve No. 

Experimental 
Graph Symbol 

Air 

Velocity 

(m/s) 

Absolute 

Humidity 

Wall 

Temperature 

(°K) 

Air 

Temperature 

(°K) 

Time 

(minutes ) 

1 

o 

1.19 

0.0070 

255 

296 

0 

1 

t 

o 

1.19 

0.0105 

255 

296 

60 

2 

A 

1.08 

0.0095 

255 

304 

0 

2 

Am* 

A 

□ 

□ 

• 

1.08 

0.0170 

255 

304 

60 

3 

1.68 

0.0093 

258 

298 

0 

3 

3.02 

0.0063 

258 

298 

105 

4 

1. 95 

0.0075 

259 

298 

0 

4 


1. 95 

0.0150 

259 

298 

60 



Theory Data 



Figure 9. Frost Thickness versus Time for Jones and Parker 
Data (Reference 4) . 


Tajima [12] for a vertical plate under laminar natural 
convection. After correcting an algebraic error (the term /Sc in 
n was placed in the denominator instead of the numerator in 
equation 23 below) in their analysis, the results are as follows, 

h H h H 

Shfi = -^ = t{l+(Us) . NUy= ^ = n/C (21), (22) 


where 


= 2/Sc 


Gr 


H 


240 [-|y + Sc$] 


1/4 


Gr = ^ 
^ H 2 

V 


■5 (T^-TJ 


f OJ — 03 1 

^ a s ^ 
[0.6453 + 2. 6453o)^J 


(23) 


(24) 


C = / + 2|^ $ + 1 - $ 

^r 

1 + OJ 

$ = ^ 

1 + 0 ) 

s 

The mean Nusselt and Sherwood numbers are 

NUjj = , Sh^ " (27), ( 28) 

The Nusselt and Sherwood numbers (Nu^ and Sh^) by these 
equations fitted Nakamura's data quite well. 

NUg and Shjj for turbulent flow, derived in Appendix III, 

are 

Sh„ = 0.029791 Sc"^/^^ 

n 


l+0.49426$Sc 


2/5 




(29) 


(25) 

(26) 


34 


c = 1.875 - /O. 765625 + 4.08973 ♦ 


(31) 


The mean Nusselt and Sherwood numbers are 


ShH = fShjj , 




(32) ,(33) 


The Olcino and Tajima [121 data show frost formations over 
time periods up to 20 hours, thus providing a severe test of the 
frost formation model. Frost was formed on vertical cylinders 
one meter in height. Frost thickness measurements were made at 
locations 0.1, 0.5, and 0.9 meters from the top. A comparison 
between equations (21) and (29) shows that the laminar flow occurs 
at 0.1 meters from the top and the turbulent flow occurs at 0.5 
and 0.9 meters from the top, given the typical input values. The 
inputs into the frost formation model are shown in Table V and 
the corresponding frost thickness data from Okino and Tajima are 
reproduced in Figure 10, The curves 1, 2, and 3 (representing 
the different measurement locations) show a very good prediction 
of the frost thickness versus time. For these curves no phase 
three regime (including water seepage) occurred. For other test 
conditions, Okino and Tajima observed cyclic frost surface tem- 
peratures indicating water seepage. When curves 4 to 9 are exa- 
mined we see there are three cycles of diffusional growth and 
water seepage growth. Curves 4, 5, 6, 7, 8, and 9 had the first 
water seepage temperatures at 268.3, 267.3, 267.3, 266.0, 264.6, 
and 265. 4“ K respectively. Thereafter the water seepage tem- 
perature was always increased by 1,2®K for successive cycles as 
corresponding to that observed by Okino and Tajima until 273.16*K 
was reached as a limit. 


TABLE V 

DATA INPUT TO THE FROST FORMATION MODEL FOR COMPARISON 
WITH OKINO AND TAJIMA [11] DATA FOR NATURAL CONVECTION ON VERTICAL CYLINDER 


Computer 
Curve No. 

Experimental 
Graph Symbol 

Plate 

Height 

(m) 

Ambient Relative 
Humidity 

Wall 

Temperature 

("K) 

Air 

Temperature 

("K) 

1 

• 

0.1 

53 

267.16 

291.16 

2 

0 

0.5 

53 

267. 16 

291.16 

3 

o 

o 

• 

VO 

53 

267.16 

291.16 

4 

A 

0.1 

64 

263.16 

292.16 

5 

A 

0.5 

64 

263.16 

292.16 

6 

A 

o 

• 

VO 

64 

263.16 

292.16 

7 

■ 

0.1 

54 

256.16 

293.16 

8 

m 

0.5 

54 

256.16 

293.16 

9 

□ 

0. 9 

54 

256.16 

293.16 


Thickness (cm) 



Time (hours) 



Frost Thickness versus Time for Okino and Tajima 
Data (Reference 11) . 


4.3 LAMINAR FORCED CONVECTION ON A CYLINDER 

Both Schneider [9] and Andrichak [13] formed frost on a 
cylinder under laminar forced convection. Since Schneider 
observed that frost thickness did not vary around the cylinder 
circumference, and since he only reported frost thickness, a mean 
Nusselt number correlation was utilized from Fand and Keswani 
[19] as. 



0.184 + 0.324Re^/^ + 0.291Re^, 


(34) 


n 1 ft — 9 ft 

X = 0.247 + 0.0407R6 for 10 < Re < 2 x 10 


The Sherwood number was obtained via equation (18). The 
data input from Schneider is shown in Table. VI and the 
corresponding frost thickness observations are shown in Figure 
11. In all cases the water seepage temperature was set at 
273..16“K. The curves 1 to 4 seem to indicate the water seepage 
model has the effect of making the frost thickness independent of 
the Reynold's number. Curve 5 predicts the data satisfactorily 
and curve 6 can be improved by using a water seepage temperature 
below 273.16°K. Schneider's correlation for the frost thickness 
(as shown by the dash lines) fits his own data quite well but 
does very poorly when applied to published data presented earlier 
in this paper. Schneider has several more frost thickness data, 
but all require water seepage temperature below freezing in order 
to obtain a proper fit with the data. Indeed, one can propose an 
empirical relationship for the water seepage temperature as a 
function of variables such as wall temperature, saturation ratio, 
etc. that would fit Schneider's frost thickness data as well as 
his correlation. The proposed empirical relationship should also 
be extended to other experimental data. 


Andrichek has a sample of frost thickness measured near the 
stagnation and near the separation points of the cylinder, as 
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TABLE VI 

DATA INPUT TO THE FROST FORMATION MODEL FOR COMPARISON 
WITH SCHNEIDER DATA [8] FOR FORCED CONVECTION ON A CYLINDER 


Computer 
Curve No. 

Experimental 
Graph Symbol 

Reynolds 

Number 

Relative 

Humidity 

Wall 

Temperature 

(“K) 

Air 

Temperature 

(“K) 

1 

0 

4000 

0. 96 

268.16 

278.16 

2 

A 

8000 

0. 96 

268.16 

278.16 

3 

0 

16000 

0.96 

268.16 

278.16 

4 

□ 

32000 

0.96 

268.16 

278.16 

5 

ID 

32000 

0. 96 

258.16 

278.16 

6 

■ 

32000 

0. 96 

248.16 

278.16 


u> 

VO 


Thickness (cm) 


Theory Data 

1 o 



Figure 11. Frost Thickness versus Time for 
Schneider's Data (Reference 8). 
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shown in Figure 12. The local Nusselt number for laminar forced 
convection on a cylinder is, 

NUj^ = 1.14Pr°'^/Re~coiT . (35) 

Equation (35) is derived from the Handbook of Heat Transfer 
[14] and fits Nusselt number data presented in Schlichting [20] 
very well up to 80° from the stagnation line. The Sherwood 
number was obtained from Equation (18) as usual. 

With the water seepage temperature set at 273.16°K it 
appears that curves 1 and 2 in Figure 12 predict the same frost 
thicknesses near stagnation and near separation. Very good 
agreement is made with the near separation data but only fair 
agreement with the near stagnation data. Yet, Schneider made 
a general observation that the frost thickness should not be a 
function of location on the cylinder. Thus the frost thickness 
distribution on the cylinder in forced convection remains 
somewhat unresolved. 

4.4 LAMINAR AND TURBULENT FORCED CONVECTION ON A FLAT PLATE 

In the last published data examined in this paper. 

White [10] used air suction through a horizontal duct with the 
frost formation on the upper duct surface. The heat and mass 
transfer coefficients were difficult to assess. For example, on 
some data sets it could not be determined whether the air flow 
was laminar or turbulent. Air suction in the duct made it dif- 
ficult to decide between treating the geometry as a duct or as a 
flat plate. It was finally decided to treat the problem as a 
flat plate, and only accept those data that show good agreement 
with frost weight predictions. The Nusselt number for the lami- 
nar flat plate is given by 

NUjj = 0.332 Rex^/^Pr^/2 (36) 
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Figure 12. Frost Thickness versus Time for Andrichak Data 
(Reference 12). (See also Parish and Sepsy) . 
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and for the turbulent flat plate given by a more complicated 
form , 

NUj{ = StjjRePr, (37) 


C /2 /UT 

-gl — = 0.9 {1 + 5/ 0. 


0897467 + ( 


Pr 

0.9 


- 1 ) 


’ 


(38) 


and 


Cf/2 = 


0.185 

(logioRe^) 


2.584 


as obtained from the Handbook of Heat Transfer. 


(39) 


The model input conditions for selected data from White are 
listed in Table VII and the corresponding frost thicknesses and 
weights are shown in Figures 13 and 14 respectively. V'Je make 
special note of curve 4 for which the frost thickness and weight 
are predicted quite well. This is the only data obtained by 
White that did not exhibit water seepage. For the first three 
curves a water seepage temperature of 27 3. 16 “K was used. Of the 
'first three curves, curve 1 provides the best prediction of frost 
thickness. Potentially, we can get good prediction of White's 
data providing we have good heat and mass transfer coefficients 
and good estimates for the water seepage temperatures below 
273.16°K. 
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TABLE VII 

DATA INPUT TO THE FROST FORMATION MODEL FOR COMPARISON 
WITH WHITE 19] DATA FOR FORCED CONVECTION ON A FLAT PLATE 


Computer 
Curve No. 

Experimental 
Graph Symbol 

Local Reynolds 
Number x 10^ 
(Laminar or Turbulent) 

Absolute 

Humidity 

Wall 

Temperature 

CK) 

Air 

Temperature 

(•K) 

1 

A 

T 

2.50 

0.00715 

242.6 

296.2 

2 

□ 

L 

1.71 

0.00762 

247.1 

296.6 

3 

A 

T 

2.68 

0.0152 

248.1 

295.9 

4 

■ 

T 

2.25 

0.0054 

250.4 

297.2 


>u 



Thickness (cm) 



Time (hours) 


Figure 13. Frost Thickness versus Time for White's Data (Reference 9) 


CSJ CO ^ 



Time (hours) 


Figure 14. Frost Weight versus Time for White's Data (Reference 9). 


SECTION 5 
CONCLUSION 

It was demonstrated in this paper that the basic approach 
by Brian et al.[l] has been improved to predict frost formation 
for other experimental data sets. Important enhancements 
included a generalized frost thermal conductivity and a water 
seepage model. Careful attention was given to deriving the heat 
and mass transfer coefficients for each particular data set. In 
particular/ for turbulent forced convection in a duct the heat 
transfer coefficient was affected by the roughness of the frost 
layer and the mass transfer coefficient was affected by the boun- 
dary layer fogging. For natural convection it was determined 
that the heat/ masS/ and momentum transfer were coupled together/ 
while the roughness effects were insignificant. 

The success of its application to a diversity of flow 
regimes/ geometrical shapes/ time dependent ambient conditions/ 
and long periods of time attest to the generality of the frost 
formation model. The frost formation model can be made even more 
general by devising an empirical relationship for the water 
seepage temperature and improving the estimates of the frost 
roughness height that can be used to calculate the heat transfer 
coefficients in turbulent flow. Eventually/ the frost formation 
model will be further modified to calculate nocturnal frost for- 
mation on a wing section. 
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APPENDIX I 

BOUNDARY LAYER FOGGING 


If the temperature in the boundary layer falls below the 
dew point, it is possible for fog droplets to be nucleated. The 
occurance of such boundary layer fogging enhances the water vapor 
mass flux to the surface. The procedure followed for calculating 
this enhanced mass flux is outlined by Epstein [21] . The mass 
transfer coefficient in the presence of fogging is given by 


.. . 1 * =p'’’a - ' =s' 

”> ” 1 + a • CpT^/L^X^ 


where 


X = 


P + P 


(40) 


(41) 


is the water vapor mass fraction and 


A = 


T. 


^v'^f 

•^a - '^f ^ (f T - I 2:3^ - 


(42) 


The fog point temperature T^ is determined by solution of the 
equation 

^ - 2.38) (T^ - Tf) ) r - 2.38) (T^ - T.)“l 

s exp |^5 ^ £-! - s. I 1 t ? ,-2 L. \ ( 43 ) 


(Lg - 2.38) (T^ - Tg) 

= 

(L^ - 2.38) (T^ - T.) 
1 + ® ^ 


1 ^ 

^ ^ 2 
L ''v'^f J 


where the critical saturation ratio for fog nucleation is 
approximated by the expression 


AnS = 2.38/R T - 2.54 
c V a 


(44) 


Equation (40) is used for the mass transfer coefficient if 

the calculated value for T_ is such that T < T^ < T . The dew- 

f s £ a 

point temperature equation is 
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where 




1 + -jp- Zn(— ) 

e ^a 


( 45 ) 


is the free stream relative humidity ratio. 
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APPENDIX II 

FROST THERMAL CONDUCTIVITY 


The frost thermal conductivity from Dietenberger [3] in 
capsule form is 

K = 1/4 f(3B - l)k- + (30 - l)k 

^ c 1 c u 

+ {[(3B^ - + (30^ - l)k„]^ + . (46) 

The upper limit conductivity expression for air bubbles and 
ice cylinders is 

k^ = (1 - B)kj^ + Bk^( upper limit) (47) 

where the thermal conductivity of air bubbles is given by; 


l=b = 


a = k . /k. 

eff-air r 


(48) 


and ice cylinders by; 

k = (1 - B)k. + Bk 


eff-air 


(49) 


Likewise, the lower limit of thermal conductivity is formed 
by an interpolation between thermal conductivities for ice 
spheres and ice planes. 

k = (1 - B)k + Bk (lower limit) (50) 

where the thermal conductivity of ice spheres is given by; 

kg = k.[3 + 2B(a - l)]/[3 - (51) 

and ice planes by; 



k . k - , 

1 eff-air 


(1 - B) k . + k.B 

eff-air i 


(52) 
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B is the proportion of the frost volume representing ice 
c 


spheres and ice planes given by, 

2 


= 13 . 6 (B 2 - Bj_)(B - Bj_) 




B - B. 


^ V®3 - ®1 - ®1. 


(B - B^) 


2(B3 - B^)(B2 - B^) 


1 + 


0.5264 


(53) 


for B > B. 


B = 0 for B < B, 
C — X 


B^ = 0.1726 (T/273.16) , 


B 2 = 0.751, 


B^ = B 2 + 0.3sin 


r^7 ^1 - (T/273.16)^ 

1 2 ' 2 
1 - (Ty273.16)"^ 


(54) 

(55) 

(56) 

(57) 


The other portion of the frost volume representing ice 
cylinders and air bubbles is given by 

©c “ 1 ” Bq . ( 58 ) 

The thermal conductivity of ice is given by 

= 630/T (59) 

and the effective thermal conductivity of air is given by, 

k ,, = k, + k„ (60) 

eff-air a v 


where 


k = 2.646 X 10 

Si 


-3 


, 1/2 


1 + 245 ,0-12/T 


(61) 


^V = T 


^s° / ^ \ / S 

\RvtV\V, 


- 1 


(62) 
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APPENDIX III 


SOLUTION OF BOUNDARY-LAYER EQUATIONS FOR FREE CONVECTION 

TURBULENT FLOW 


The integrated boundary-layer equations reproduced from 
Nakamura [18] for free convection are, 

(a) Momentum Equation 


d 5 2 ^ 

g- [ / pu dy] = - + / (p - p^)g cos<).dy 

o o 


(63) 


(b) Energy Equation 


puOp(T - T^jdy] = Jq„ + »Vp<'^w - ■'-> 


(64) 


(c) Diffusion Equation 


dx 


[ / pu(X^ - X^.)dy] = + pV„(X^„ - Xj^.) 


(65) 


(d) Normal velocity at the frost surface (zero net flow 
of air) 


'Iw 


" p(l - X^;^) 


( 66 ) 


Barron and Han [22] show the same equations with one 
important subtle difference. Barron and Han assumed the thermal 
thickness is the same as the momentum thickness, 6, and the dif- 
fusional thickness is different from the momentum thickness. 
Nakamura [18] and Skelland [23] assumed that the diffusional 
thickness is the same as the momentum thickness and the thermal 
thickness is different from the momentum thickness. We chose the 
latter approach for the reason explained later. Barron and Han 
demonstrated the effect of thermal diffusion as nil, thus 
simplifying the equations somewhat. The integrals were evaluated 
using the following velocity, temperature and mass fraction pro- 
files from Skelland and Barron and Han as 
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u. 




m 


0 = 


where 


U = 
m 

6 = 
C = 


P oo 


®t T 


and 


B = 


0.6078 


( 67 ) 


T - T“ 

1/7 

= 2 - . 

(68) 

T — T» 
w 

Xi - Xj_ 

. = [1 - (/l„ - Xl„) 

(69) 

C^X^/2 

f 

(70) 

c.xVio 

r 

(71) 

St/'S / 

and 

(72) 

-B^(T - 

T«) - fi^(X^ - X^„) 

(73) 


m 1 + 0.6078X. 


(74) , (75) 


From Skelland the surface shear stress was 

u 2 1/4 

T = 0 . 0228 p(-|-)(^) 

^ c m 


(76) 


Then the Colburn's analogy was used to get the heat and mass 
transfer equations as , ■ 


^Iw/Vc == 


^iw - ^1.. 

m 


and 


1 /T g = 
-'qw' o^c 


Pr^/^U 


(77) 


(78) 


m 


Substituting all the previous equations into equations 
(63), (64), and (65) the following result is obtained. 
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088935 = S_£|S<|. [ 58 ^( 1 ^ - T^) + - X^^)] 


-0.0228V 


1/4 


7\ 1/4 


( 79 ) 


- 1/4- 5/4 6 r.8/7,7 7. .15/7.^ - 28 W ^22/7 .42 _ 42 

'“2 5 *^8 9-< ^ '•15 le-* ^ '•22 23 ^ 

29/7.^ . p36/7.7_ _ 7_^, 

" ^ '•29 30-' ^ '•36 37 


^ ^ — 173 ■'■ 

Pr 


X, - X, 
Iw 1< 


Sc2/3(i - 


) [ 0 . 0228 v^/^) 


(80) 


|(0. 03663) 

- {[ 1 )[0.0228v)^/^ (81) 

The solution for 5 is obtained from combining equations 
(80) and (81) and doing a parabolic Taylor's expansion of the 
terms containing about C = 1 to get equation (31) in the text. 

The values for and C 2 are solved from combining 
equations (79) and (80) which then lead us to the equations, for 
6 and U as. 


X" 


1 “ X, -0.8 - , - 


1 - X, ,,, 0.1 
(1 + 0.4943( j - ■ - - j^- -)Sc^^'^) 


(82) 


1-X 2/3 —1/2 

and = 1.1855(.^) 0.4943 )Sc ) 

loo 


(83) 


where Gr^ is given by equation (24) in the text and 
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I 


. - ^ 1 * “a 

V I XT ”1 i - . f 


1 - X, 


1 + U). 


1 - X, 


= 1 + 0)^ , H = 


•* "s "Iw 

The Sherwood and the Nusselt numbers are defined 




and Nu^ - 


X. (84) 
by 

(85) ,(86) 


which when combined with equations (76), (77), (78), (82), and 
(83) give the equations (29) and (30) in the text for the local 
Sherwood and Nusselt numbers. As a result we obtained the 
Chilton-Colburn analogy given by equation (30) as should happen 
when using Colburn's equations for the heat and mass transfer. 
Barron and Han derived a different analogy result than the 
Chilton-Colburn analogy, thus invalidating their own results. 

To evaluate the roughness effects on the heat transfer 

coefficient in turbulent flow the friction velocity, U^, obtained 

from the stress equation (76) at a one-meter height for typical 

conditions was combined with the typical frost thickness to show 

that the roughness Reynolds number, U k /v was around 5 to 25, 

1 s 

thus indicating a smooth wall behavior for the Okino and Tajima 
data. 
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